Counterion Condensation and Fluctuation-Induced Attraction 
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We consider an overall neutral system consisting of two similarly charged plates and their op- 
positely charged counterions and analyze the electrostatic interaction between the two surfaces 
beyond the mean-field Poisson-Boltzmann approximation. Our physical picture is based on the 
fluctuation-driven counterion condensation model, in which a fraction of the counterions is allowed 
to "condense" onto the charged plates. In addition, an expression for the pressure is derived, which 
includes fluctuation contributions of the whole system. We find that for sufficiently high surface 
charges, the distance at which the attraction, arising from charge fluctuations, starts to dominate 
can be large compared to the Gouy-Chapmann length. We also demonstrate that depending on the 
valency, the system may exhibit a novel first-order binding transition at short distances. 

PACS numbers: 82.70.-y,61.20.Qg 
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I. INTRODUCTION 

Correlation effects may play an important role in 
controlling the structure and phase behavior of highly 
charged macroions in solutions The macroions may 
be charged membranes, stiff polyelectrolytes such as 
DNA, or charged colloidal particles. Recently, these ef- 
fects have attracted a great deal of attention, since they 
may drastically alter the standard mean-field Poisson- 
Boltzmann (PB) picture [|, U |, |, §. For example, 
one surprising phenomenon is the attraction between two 
highly-charged macroions, as observed in experiments 
@, I I and in simulations [|o[ This attrac- 

tion is not contained in the mean-field (PB) theory, even 
for an idealized system of two charged planar surfaces. 
Indeed, it has been proven that PB theory predicts only 
repulsion between likely-charged macroions jT^ . 

Very recently, another interesting effect that is not cap- 
tured within the PB theory is predicted, namely the novel 
fluctuation-driven counterion condensation [Q. For a 
system consisting of a single charged surface and its op- 
positely charged counterions, Netz and Orland Q showed 
that a simple perturbative expansion about the mean- 
field PB solution breaks down for sufficiently high sur- 
face charge. Thus, in this limit, fluctuation and cor- 
relation corrections can become so large that the solu- 
tion to the PB equation is no longer a good approxi- 
mation. To circumvent this difficulty, a two-fluid model 
was proposed in Ref. [l4| , in which the counterions are 
divided into a free and a condensed fraction. The free 
counterions have the usual 3D mean field spatial distri- 
bution, while the condensed counterions are confined to 
move only on the charged surface and thus effectively re- 
duce its surface charge density. The number of condensed 
counterions is determined self-consistently, by minimiz- 
ing the total free energy which includes fluctuation con- 
tributions. This theory predicts that if surface charge 
density of the plate is sufficiently high, a large fraction of 
counterions is "condensed" via a phase transition, simi- 
lar to the liquid-gas transition with a line of first-order 



phase transitions terminating at the critical point. Fur- 
thermore, the valence of the counterions plays a crucial 
role in determining the nature of the condensation transi- 
tion. The physical mechanism leading to this counterion 
condensation is the additional binding arising from 2D 
charge-fluctuations, which dominate the system at high 
surface charge. In this paper, we extend this conden- 
sation picture to a system of two charged surfaces with 
their neutralizing counterions, and to study the electro- 
static interaction between them. 

Previous theoretical approaches to the problem of the 
attraction in charged surfaces includes both numerical 
and analytical methods that go beyond the mean-field PB 
theory. Gulbrand et al. provides the first convincing 
demonstration for the attraction between highly-charged 
walls using Monte Carlo simulations. In particular, they 
showed that for divalent counterions, the pressure be- 
tween charged walls becomes negative for distances less 
than 10 A; Hence, the existence of short-ranged attrac- 
tion. Subsequently, there has been a number of numerical 
studies based on the hypernetted chain integral equation 
JlSf and the local density- functional theory |l^, as well 
as analytic perturbative expansion around the PB solu- 
tion that demonstrate attraction. 

More recently, motivated by the problems of DNA 
condensation and membrane adhesion, two distinct ap- 
proaches have been proposed to account for the attrac- 
tion arising from correlations |^, |4j. The first ap- 
proach based on "structural" correlations first proposed 
by Rouzina and Bloomfield ||^ , the attraction comes from 
the ground state configuration of the "condensed" coun- 
terions. This theory predicts a strong short-ranged at- 
traction, with the characteristic length set by the lattice 
constant, typically of the order of few Angstroms. In the 
other approach, based on charge-fluctuations, the coun- 
terion fluctuations are approximated by the 2D Debye- 
Hiickel theory, which predicts a long-ranged attraction 
which scales with the interplanar distance as d~^ . Note, 
however, that the mean-field PB repulsion which scales 
like d~'^ always dominates the attraction for large dis- 
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tances, and thus, the range of the attraction in this pic- 
ture is still short, typically of the order of 10 A|l^, p^ . 
Despite of the fact that some conceptual issues have been 
resolved concerning the crossover of the attractions from 
long-ranged to short-ranged |2^, there remain some 
interesting problems to be understood. In particular, 
some experimental observations in planar surfaces |^ and 
in charged colloidal suspensions |9[] as well as computer 
simulations provide evidence for a long-ranged at- 
traction, typically of the order of microns, whereas the 
two mechanisms mentioned above give only short-ranged 
attraction. In this paper, we show that the charge- 
fluctuation approach, together with the counterion con- 
densation mechanism , can induce long-ranged attrac- 
tions for sufficiently high surface charge. We note that 
other mechanisms based on hydrodynamic interactions 
p2|, depletion effects 1231, and an exact calculation for 



the 2D plasma model |24 have been proposed recently 
to account for the long-ranged attractions. 

In particular, we study the interaction between two 
charged surfaces separated by a distance d, with counte- 
rions distributed both inside and outside of the gap. This 
boundary condition, as opposed to all of the counterions 
confined between the gap, is more appropriate in general, 
since systems are not closed and often the counterions 
are in equilibrium with a "bath" in surface forces experi- 
ments. In the spirit of the "two-fluid" model proposed in 
Ref. we divide the counterions into a "condensed" 

and a "free" fraction. The condensed counterions are al- 
lowed to move only on the charged surfaces, while the free 
counterions distribute in the space inside and outside the 
gap. The surface density of the condensed counterion Uc 
on each plate is determined by minimizing the total free 
energy, which includes fluctuation contributions. Fur- 
thermore, an expression for the fluctuation pressure is de- 
rived, which includes fluctuation contributions from the 
condensed and "free" counterions, and their couplings. 
We find that the counterion condensation can occur ei- 
ther by increasing surface charge density at a fixed dis- 
tance or by decreasing the separation between plates. For 
low surface charge, the counterion condensation proceeds 
continuously as a function of distance with the fraction 
of counterion condensed being small but finite, and the 
total pressure of the system remains repulsive. 

For higher surface charges, the qualitative behavior of 
the counterion condensation transition depends critically 
on the valence Z of the counterions. For Z < 2, the coun- 
terion condensation proceeds continuously as a function 
of distance. On the other hand, for Z > 2, the behav- 
ior of the system is qualitatively different, similar to an 
isolated charged plate |Q. In this case, the counterion 
condensation occurs via a first-order phase transition as 
a function of distance. Remarkably, we find that for 
trivalent (Z — 3) counterions, there is a wide range in 
the surface density, in which the first-order counterion 
condensation spontaneously takes that the system from 
a repulsive regime to an attractive regime at short dis- 
tances, resulting in a first-order binding transition. For 



high surface charge, counterion condensation again pro- 
ceeds continuously even for Z > 2, but with a significant 
number of condensed counterions. Thus, in this regime, 
the mean-field repulsion is substantially reduced and the 
long-ranged charge-fluctuation attraction dominates the 
system even for large distances. Note, however, that the 
mean- field repulsion will eventually dominate as d — > oo. 
We emphasize that all these features, in particular the 
special role of the valence, deviate significantly from the 
PB mean-field predictions. 

This paper is organized as follows: In Sec. we briefly 
recapitulate qualitatively the mechanism which drives 
the counterion condensation. In Sec. |l|, we present 
in detail the two-fluid model and derive a general expres- 
sion for the total free energy of the counterions. In Sec. 
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we apply this formalism to study the interaction of 
similarly charged surfaces. A detailed discussion of our 
results is presented in Sec. 0. 



II. COUNTERION CONDENSATION: 
QUALITATIVE ARGUMENT 

In this section, we recapitulate the essential physics of 
the condensation transition presented in Ref. [Q . Recall 
that for a single plate of charge density cr(x) = eno(5(z) 
immersed in an aqueous solution of dielectric constant 
e, containing oppositely charged —Ze point-like counte- 
rions of valence Z on both sides of the plate, PB theory 
predicts that the counterion distribution psl 
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decays to zero algebraically with a characteristic length 
A = I/^ttIbZhq), where Ib = j^-^ ~ 7 A is the Bjerrum 
length in water at room temperature, fc^ is the Boltz- 
mann constant, and T is the temperature. This Gouy- 
Chapman length A defines a sheath near the charged sur- 
face within which most of the counterions are confined. 
Typically, it is on the order of few angstroms for a mod- 
erate charge density of no ^ 1/100 A~^. Note that since 
A scales inversely with hq and linearly with T, at suffi- 
ciently high density or low temperature, the counterion 
distribution is essentially two-dimensional. In fact, in the 
limit T 0, we have 



lim / pn(z)dz — lim 2 • — — 
r^Q ^ T^o 2Z 



C 



A 
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where C is an arbitrarily small but fixed positive value 
of z, i.e. the counterion profile pq{z) reduces to a sur- 
face density coating the charged plane with a density of 
hq/Z. Therefore, according to PB theory, all of the coun- 
terions collapse onto the charged plane at zero temper- 
ature. However, for highly charged surfaces Z'^Ib ^ A, 
the fluctuation corrections become so large that the solu- 
tion to the PB equation is no longer valid . To capture 
this regime in the spirit of the "two-fluid" model p4|. 
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FIG. 1: The geometry of the problem. 



we divide the counterions into a "free" and a conden- 
sate fraction. The "free" counterions have the usual PB 
3D spatial distribution, while the "condensed" counteri- 
ons are confined to move only on the charged plane, as 
shown in Fig. ||. The free energy per unit area for the 
condensed counterions with an average surface density ric 
can be written as [Eql 



fihd{nc) 
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where (3~^ — ksT, a is the molecular size of the counte- 
rions, and Xd — l/(27rZ^nc) is the 2D screening length. 
The first term in Eq. is the entropy and the sec- 
ond term arises from the 2D fluctuations. Note that 
the latter term is logarithmically divergent, which may 
be regularized by a microscopic cut-off a, yielding 
PAf2d(nc) — ~ gJx'^ ln(27rAjj/Q). In addition, the con- 
densate partially neutralizes the charged plane, effec- 
tively reducing its surface charge density from eng to 
enji = eriQ — Zen^. Thus, motivated by PB theory, the 
free counterions can be modelled as an ideal gas confined 
to a slab of thickness A^, = Xj [t^IbZub) with a 3D con- 
centration of c = n^ jiZXji). The fluctuation free energy 
in this case may be estimated using the 3D Debye-Hiickel 
theory: /3A/ = -K3/(127r) ||||, where kI = A.txZHbc is 
the inverse square of the 3D screening length. The free 
energy per unit area of the free counterions is then ap- 
proximately given by 



(4) 



All the qualitative results, including the nature of the 
condensation transition, follow straightforwardly from 
the analysis of the total free energy: /(ric) = /2d("-c) + 
fadijic)', minimizing f{nc) to find the fraction of con- 
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where the three dimensionless parameters: the order pa- 
rameter T = Zric/no, the coupling constant g = Z^Ib/X, 
(where A is the hare Gouy-Chapmann length), and the 
reduced temperature 9 = ; completely determine 
the equilibrium state of the system. It is easy to de- 
rive the asymptotic solutions of the last equation corre- 
sponding to the free, ri <C 1, and condensed, T2 ~ 1, 
state of the counterions: ri — gO exp [l — 1 5] j and 

g-l 

T2 — 1 — [7rexp(l)] ^ , respectively. For weak 

couplings g ^ 1, Ti is the only consistent solution. Thus, 
there is almost no condensed counterions r ^ 1. This is 
not surprising since PB theory is a weak-coupling the- 
ory which becomes exact in the limit T ^ 00. However, 
for high surface charge g S> 1, where correlation effects 
becomes important, the behavior of t depends crucially 
on 9. In particular, for 9 < 9c ^ 0.038, ti and T2 are 
both consistent solutions corresponding to the two min- 
ima of /, and thus a first-order transition takes place 
when /(ti) = /(''2), in which a large fraction of con- 
terions is condensed. This occurs at a particular value 
of the bare surface charge density such that g = go{9). 
For an estimate, we take = 0.02 (divalent counterions 
at room temperature) and obtain go 1.7, correspond- 
ing to CTc ^ e/10 nm^^. However, for 9 > 9c the be- 
havior of T is completely different; in this regime, there 
is no phase transition and the condensation occurs con- 
tinuously. Thus, the condensation transition is similar 
to the liquid-gas transition, which has a line of first- 
order transitions terminating at the critical point where 
a second-order transition occurs. If one takes Ib ^ W 
A, e.g. room temperature, and a ~ lA, it follows from 
the definition of 9 that there is a critical value of coun- 
tcrion valence Zc — \/a/{lB9c) — 1.62, below which no 
first-order condensation transition is possible. Therefore, 
divalent counterions behave qualitatively differently from 
monovalent counterions at room temperature. 

Clearly, this condensation picture may also be cru- 
cial to understanding the attraction between two simi- 
larly charged plates, separated by a distance d. Recall 
that the total pressure of this system is comprised of the 
mean-field repulsion and the correlated fluctuation at- 
traction Q. The repulsion comes solely from the ideal 
gas entropy and it is proportional to the concentration 
at the mid-plane: no(d) = ksTpoiO) = 8/cbT/(£bA|j) 
for d < Xii The fluctuation-induced attraction is 

n(d) = -aokBT/d^ for d > Xd, where ao « 0.048 
Q. Clearly, when a large fraction of the counterions is 
"condensed" , the mean-field repulsion is greatly reduced. 
Therefore, the attraction arising from correlated fluctu- 
ations of the "condensed" counterions can overcome the 
mean-field repulsion even for large distances. Using the 
estimates in the last paragraph above, we find that for 
divalent counterions and surface charge density of about 
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one unit charge per S ~ 7nm^, the total pressure be- 
comes attractive at about d 10 nm; hence a long- 
ranged attraction. Of course, this estimate should be 
supplemented by a more precise calculation for the sys- 
tem of two charged plates, which is done below. 

III. COUNTERION FREE ENERGY IN THE 
"TWO-FLUID" MODEL 

Consider an overall neutral system consisting of coun- 
terions and two charged surfaces separated by a distance 
d immersed in an aqueous solution. The surface charged 
density on each plate is (Tq = eriQ- We model the aque- 



ous solution with a uniform dielectric constant e. This 
simplification allows us to study fluctuation and correla- 
tion effects analytically. In the spirit of the "two-fluid" 
model proposed in Ref. IT^ , we divide the countcrions 
into a "condensed" and a "free" fraction. The condensed 
counterions are allowed to move only on the charged sur- 
faces, while the free counterions distribute in the space 
inside and outside the gap. The number of the condensed 
counterion ric on each plate is determined by minimizing 
the total free energy including fluctuation contributions. 
Thus, our first task is to derive an expression for the total 
free energy of the system. The electrostatic free energy 
for the whole system may be written as 



l3Fei = J2 {In K(r)a'] - l} + j d^^p{^) {in [p{^)a^] - l} 



ZIbY^ J d^x J d^x' 



,3 , nl{r)6{z - Zj) nl{v')6{z - zj) Z'^li 



|x-x'| 

,n^(r)<5(2-z,)[Mx')-"/(x')] 



d^x / d^xl 



, P(x)p(x') 



(6) 



where a is the molecular size of the counterions, Ib = 
e^/ (eksT) is the Bjerrum length, Z is the valence of the 
countcrions, and Zi = —d/2 and Z2 = d/2 are the loca- 
tions of the charged surfaces. The first two terms in Eq. 

are the two-dimensional entropy for the condensate 
and three-dimensional entropy for the "free" counterions, 
respectively, and the rest represent the electrostatic in- 
teractions of counterions in the system. In Eq. (^, the 
condensed counterions two-dimensional density on the i- 
th plate is denoted by nKr), the "free" counterions with 
3D density by p{x), and the external fixed charges arising 
from the surfaces by (x) = no6{z — d/2)+noS{z + d/2). 
Within the Gaussian fluctuation approximation, we as- 
sume that the 2D density of condensed counterions has 
a spatially dependent fluctuation about a uniform mean: 
(r) = 7ic + Snl{r), and expand Eq. (||) to second order 
in 5ri*(r). Summing over all the 2D fluctuations of the 
condensed counterions, i.e. 



-/3-He 



we obtain two terms in the effective free energy: Tie = 
F2d+Ti-j,d- The first term F2d is the free energy associated 
with the condensed countcrions which can be written as 



2d 



2nc {in 



1 



ilndct[-V^], 



(7) 



where K2rf(x,y) = [-V^ + 5(z ± d/2) J 5(x - 

y) is the 2D Debye-Hiickel operator and \d = 
\/{2TTZHBnc) is the Debye screening length in 2-D. The 
first term in Eq. is the entropy and the second term 
arises from the 2D charge-fluctuations. Note that al- 
though this fluctuation term can be evaluated analyti- 
cally [Q, we write it in this abstract form for later con- 
venience. 

The second term Ti^d is the electrostatic free energy for 
the "free" counterions, taking into account of the pres- 
ence of the fluctuating condensate; to within an additive 
constant, it may be written as 

m-M ^ y"rf'xp(x){ln[p(x)a3] -1} 

+ ^ Jd^xJ dVp(x)G2d(x,x')p(x') 

d^x(t,{x)p{x), (8) 



Tic a 



1}^ + - IndetKzd 



where (j){x) = J d^x! Z^^ G2ci(x, x') nfl(x') is the "renor- 
malized" external field arising from the charged plate. 
From Eq. (^) , we can see that the presence of the conden- 
sate modifies the electrostatics of the free counterions in 
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two ways. First, the condensate partially neutralizes the 
charged surfaces, effectively reducing the surface charge 
density from eng to enji = e(no ~ Zuc). Second, their 
fluctuations rcnormalizc the electrostatic interaction of 
the system; thus, instead of the usual Coulomb poten- 
tial, the free counterions and the charged surfaces inter- 
act via the interaction G2(j(x, x'), which is the inverse 
(the Green's function) of the 2D Debye-Hiickel operator 
K2d: 



+ — ^ 6{z ± d/2) G2<i(x, x') = £bS{^ ~ x'), 

(9) 

where we have defined, for convenience, a reduced Bjer- 
rmii length hy is = A-kZ^Ib- In Eq. (||), the sec- 
ond term in the bracket takes into account of the fluc- 
tuating "condensate" . Hence, in the limit Uc or 
Xd oo, G2d(x, x') reduces to the usual Coulomb inter- 
action Go(x, x') — £b/\^ — x'|. 

After a Hubbard-Stratonovich transformation the 
grand canonical partition function for the "free" counte- 
rions, characterized by the interaction energy in Eq. (|^), 
can be mapped onto a functional integral representation: 
Zf,[(t>] = A/'o / X>i/;e-'5['^>"^l with an action 



- J ^ |ivXx)[-V^]V.(x)-.^e 
+ -L^<5(z±d/2)[V;(x)]4 



2 j^(x)-|-0(> 



(10) 



where tp{x.) is the fluctuating field, ~ e^^s/a^, fi is 
the chemical potential and A/q"^ = det 'K2d is the nor- 



malization factor. The minimum of the action, given by 
-g^^ — 0, defines the saddle-point equation for 

■00 (x) 



-V2[iVo(x)] + — V<5(z±d/2)[z^o(x)] 

-f At2g#o(x)+0(x) ^ Q 

This saddle-point equation is equivalent to the PB 
equation by defining the mean-field potential ^(x) — 
— zi/'o(x) — 0(x), which is solved below in Sec. IV. To 
obtain the free energy for the free counterions on the 
mean-field level, we note that it is related to the Gibbs 
potential To[(j)\ = SlipOi by a Legendre transformation: 



PF^^dinn) = ToM -f d3xpo(x). 



(12) 



where po(x) is the mean-field free counterion density 
given by 



po(x) = {h^/£b) e*'''°(")+"*(") . 



(13) 



To capture correlation effects, we must also include 
fiuctuations of the "free" counterions, thereby treating 
the "free" and "condensed" counterions on the same 
level. To this end, we expand the action S[ip, 0] about 
the saddle-point '!/'o(x) to second order in A'0(x) = 
V'(x) - '00 (x): 



S[^, 0] = 5[0, ^o] + 2 / I d y A0(x) K3<i(x, y) Ai^{y) + ■■■, 



where the differential operator 
K3d(x,y) 



^(5(z±d/2) + K 

± 

I 



2giVo(x)+0(x) 



-^(x-y), 



(14) 



is the second variation of the action S[ip,(t)]. Note that 
the linear term in Aip{yi) does not contribute to the ex- 
pansion since 0o(x) satisfies the saddle-point equation 
Eq. ([Tl] ) . Performing the Gaussian integrals in the func- 
tional integral, we obtain an expression for the change in 
the free energy due to fiuctuations of the free counterions: 



where the second term comes from the normalization fac- 
tor Aq. Note that the second term in Eq. ( p^ ) partially 
cancels the fiuctuation contributions to the free energy of 
the condensed counterions in Eq. (Q). Thus, combining 
Eqs. (0), ([l^), and ( p^ ) together, the total free energy of 
the system can be expressed as 



PAFsd = ^ In det K^d - ^ In det K2d, (15) 



1 



-f - In det Ksd 



■lndet[-Vi]. 



(16) 
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This is the main result of this paper, from which all the 
equilibrium quantities can be calculated. It says that 
the free energy of the counterions is simply a sum of 
the mean-field free energy and a fluctuation energy term. 
Note that the latter term contains couplings among the 
fluctuations of the free and the condensed counterions. 
Finally, we stress that the derivation presented here is 
rather general and may apply to other physical systems 
as well. 



IV. INTERACTION BETWEEN TWO 
SIMILARLY CHARGED SURFACES 

In this section, we employ the framework of counte- 
rion condensation derived in Sec. 



HI to study the inter- 



action of two charged surfaces with free counterions and 
condensed counterions fluctuating on each of them. In 
Sec. 



IV A 



an expression is derived for the total pres- 
sure, which takes into account the total fluctuations of 
the counterions. 



In Sec. IV B, wc discuss the behavior of 



the total pressure and the equilibrium state of the system 
as characterized by the fraction of condensed counterion 
r which is determined by the minimum of the total free 
energy Eq. (|l|). 



A. Mean-Field Theory and Fluctuation Corrections 

The free counterion density on the mean-field level can 
be obtained by solving Eq. Defining the mean- field 

potential by ip{x) — —i^o{x) — 0(x), the saddle-point 
equation becomes 



— J2Siz±d/2Mz), (17) 



where enn = e(no — Zuc) is the renormalized surface 
charge density of the plates. Note that Eq. ( [TtI ) looks 
similar to the mean-field PB equation. Indeed, the so- 
lution to Eq. ( [l7| ) is exactly the same as the PB solu- 
tion provided that we impose the boundary condition: 
</j(±(i/2) = 0. The solution reads 



'P<{z) 



, K^cos^(az) , , 
21n l + -^(|z|-rf/2) 



(18) 

\z\ > d/2, (19) 



and the counterion density is given by 



poe = {2a^/lB)sec^{az),\z\ <d/2, (20) 



i\z\-d/2 + a' 



\z\ > d/2, (21) 



where ^ = ^/2/k and a is determined from the boundary 
conditions on the electric field: dzip>\j^^2 ~ ^^'P<\d/2 ~ 
riRiB/Z and ip{± d/2) — 0; they lead to a transcendental 
equation for a: 



uXr tan {ad/2) 1 - {aXii/2f 



(22) 



where A/j AZ / {isnR) is the renormalized Gouy- 
Chapmann length. Physically, is proportional to the 
free counterion density at the mid-plane Pq{0). In ad- 
dition, K is related to the zeros of the potential, i.e. 
ip{±d/2) = 0: 



= 2a^ scc^ {ad/2) = 2(1 + 6^) 



(23) 



where we have defined h = a\if/2. The asymptotic be- 
haviors for b as determined by the relation Eq. (|2|) are 
h ^ 1/dasd^oo and 6 ~ 1 as d ^ 0. Wc note that the 
latter behavior is distinct from the case of two impene- 
trable charged hard walls |2^ . 

The mean-field free energy per unit area for the free 
counterions, i.e., the first two terms in Eq. ([l^), can be 
easily calculated by using Eq. (12), 



Ph - 2n,{lnKa2]_i} + ^| 
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1 + iaXu/2)' 



2a^d 
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where we have made used of the fact that the chemi- 
cal potential is given hy fi = \n {k^g^ /£b) ■ The first 
two terms in Eq. (24) represent the entropy per unit 
area of condensed and free counterions, respectively. The 
last two terms describe the interacting free energy for 
the surfaces. Using the general formula for the pressure: 



dfojd) 

tween the surfaces: IlQ{d) = -^2c?/Ib. We note that 



no(d) = — "^Q^^^ , we obtain the mean-field pressure be- 



at the mean-field level, the pressure comes solely from 
the ideal gas entropy of the "free" counterions, and it 
is proportional to their density at the midplane, a stan- 
dard result. However, in contrast to the standard PB 
theory, the pressure now depends on the order parame- 
ter r = Zric/riQ. Thus, if there were a large fraction of 
condensed counterions, r ~ 1, the mean- field repulsion 
would be drastically reduced. 

Next, we compute the pressure arising from the coun- 
terion fluctuations. Recall that the expression for the 
change in the free energy arising from fluctuations of the 
counterions is given by the last two terms in Eq. ( p^ ) 



(3 AT = i In det K^a ^ ^ 1^ det 



(25) 



where the operators K-^d is defined in Eq. ([14) with 



^2g-y(2) =2a^sec^{az)e{z) + 



2e(z) 



i\z\-d/2 + a' 



(26) 



where ^ = Xr/{1+P), e{z) = 0{z+d/2)-e{z-d/2) = 1, 
if |z| < d/2 and zero, otherwise, and 0(z) = 1 — Q{z). 
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The derivative of AjF with respect to distance d can be 
straightforwardly calculated by making use of the exact 
identity: S In det ± = Tt±~H±, 



dd 



^ y"d3xG3rf(x,x) X 



d_ 

dd 



A 



D 



S{z±d/2)+ K^e 



-viz) 



(27) 



where G3rf(x, x') is the Green's function of the operators 
satisfying 



^(5(2±d/2)+ K 



2g-v(2) 



G3d(x,x') 



^B^(x-x'). (28) 



An explicit derivation of G^dhi., x') and the pressure aris- 
ing from fluctuations Eq. ( p7| ) are detailed in Appendix 
The final result for the pressure can be written as 



A dd 

(f + 



(2^)2 

I3) 



qM'^iq) 
l-M^{q) 



(29) 



where M{q) is defined in Eq. ( A14 ) , X2 and X3 are two di- 
mensionless integrals defined in the Appendix ^ by Eqs. 
( A3C| ) and ( A31 ), respectively. 

We note that the fluctuation pressure is purely attrac- 
tive; thus, fluctuations lower the free energy. Although 
the expression Eq. (^) looks complicated, each term, 
however, has a simple physical interpretation. The first 
term in Eq. (|2^) is the pressure arising from countcrion 
fluctuations near the surfaces. In fact, if all of the coun- 
terions is condensed t = 1, we observe that M{q) in Eq. 
( [A14| ) becomes M{q) = -e"«''/(l + gAu^and that the 
only contribution to the pressure in Eq. (E3) is the first 
term, which becomes in this limit 



Pll{d) 



d^q q 
(27r)2 e^i'^il+qXay - I' 



(30) 



This expression is exactly the pressure derived in Ref. 

arising from 2D fluctuations of the counterions. It 
scales like Ii{d) ^ —1/d^ for large distances. 

The second term in Eq. ( [29| ) may be interpreted as the 
coupling between the counterions near the surfaces and 
those in the bulk. This can be seen by considering the 
asymptotic behavior of the pressure for large d in the no 
condensate limit t = 0, i.e. the fluctuation corrections 
to the PB pressure. In this case, in addition to the usual 
d^^ scaling law arising from counterion fluctuations near 
the surfaces, the second term in Eq. ( |2^ ) contributes a 
term, which scales as ~ d~^ In (d/A) in the large d limit. 
Therefore, the pressure 



n(d) 



1 

1? 



d3 



In (d/A) 



(31) 



contains a logarithmic term, which dominates the d~^ 
term for large distances. This term has been obtained 
by several authors previously ||l^, |l^ and in particular, 
Ref. [|l8| shows that this term arises physically from the 
coupling between counterions near the surfaces and those 
in the bulk. Therefore, Eq. ([2S| ) recovers the PB limit 
T — Q and the 2D limit r = 1 as special cases. Although 
the fluctuation corrections to the PB (r = 0) pressure 
have been considered previously, we stress that Eq. ( [29| ) 
is a generalization which allows for counterion condensa- 
tion and may apply to other physical situations, such as 
ions absorption. 

Combining with the mean-field pressure, we obtain the 
total pressure 



mot(d) 



^ 1 ^ 

Ib \ SttXr 2 + ^(1 

d^q qM'^iq) 
(2^ 1-MHq)' 



62) 



(32) 



The behavior of the total pressure depends on the cou- 
pling constant g = Z^Ib/X and the fraction of condensed 
counterions r = Zric/nQ. For g ^ 1 and r <C 1, the 
fluctuation corrections are small and the total pressure 
^tot{d) is controlled by the mean- field repulsion. On the 
other hand, for r ~ 1 the mean-field repulsion is greatly 
reduced and the fluctuation attraction can overcome the 
repulsion at finite distances. Furthermore, for g ~ 1 the 
short distance behavior is highly sensitive to r: Even a 
very small number of condensed counterions would turn 
the total pressure, otherwise repulsive for t = 0, into 
attractive for short distances. For g ^ 1, the fluctua- 
tion attraction becomes dominant at short distances even 
when there is no condensate, and the effect of finite t is 
to push the attractive region out to a larger length scale. 
Hence, if there is sufficient number of condensed counte- 
rions, the pressure is attractive even for large distances. 
Our next task is to determine the fraction of condensed 
counterions t as determined by the minimum of the total 
free energy. 



B. Equilibrium Properties 

The equilibrium state of the system is determined by 
minimizing the total free energy with respect to the order 
parameter t. Therefore, we need to evaluate the deriva- 
tive of the total free energy Eq. (|l^ ) with respect to r. 
Let us first consider the mean-field contribution. Explic- 
itly differentiating Eq. (kJ), we obtain 



mo 

dT 



^ln(2A/a)-^ + ^ln^-^. 
Z ^ ' ' Z Z (l-r)2 

^ In [1 + 6^] , (33) 



where A = 4Z/(£b?^o) is the "bare" Gouy-Chapman 
length. To obtain the fluctuation contributions, we again 
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make use of the exact identity: 5 IndetX = TrX ^ (5X 
to evaluate the derivative of the fluctuation free energy 
in Eq. (H), 



2ii 
d_ 



^ J d^xG3<j(x,x) X 



•(34) 



This expression can be exphcitly evaluated using similar 
techniques outlined in Appendix]^ and the result is given 

by 



1 dpAT _ 4(1 - r) 
A dr " A2 



2&2(X2-X3) 



2+^(1 + 62) 



, (35) 



where T2 and X3 are given in Eqs. 
respectively, and Ii is defined by 



(A30) and (A31), 



(36) 



where £(1 7) and Q{d/2) are defined by Eqs. ( [Al8| ) 
and ( [A19| ), respectively, in Appendix Note that 
Ji[d/ Afl] is logarithmically divergent (see Appendix 0, 
Eq. (A34)), as in the case for 2D Debye-Huckel theory, 
which may be regularized by a microscopic cut-off, cho- 
sen to be the size of the counterion a. Finally, using Eqs. 
( ^ ) and (^5|), the root of the free energy dF{T)/dT = 
can be determined numerically. For example, the case of 
an isolated charged plate can be obtained by taking the 
limit d ^ 00 in Eqs. ( ^3| ) and (|3^), which leads to the 
following transcendental equation 



1 + ln(.g6l/2) -l-ln 



X dx 



l + 2-f{l + x) 



0,(37) 



where Xc 



{l+x)[l + {-f + x){l + x)^ 
is the microscopic cut-off. The anal- 



ysis of this equation gives all the features mentioned in 
Sec. ||. As a consistency check, it can be verified that 
in the limit d — *■ 0, Eq. ( p5[ ) gives the fluctuation free 
energy for an isolated charged surface but with twice of 
the surface charge density 2enQ. 



V. RESULTS AND DISCUSSION 

Let us first discuss the behavior of the order parameter 
r at a fixed separation d between the charged surfaces, 
as summarized in Fig. |^. The behavior of r as a func- 
tion (7 at a finite distance d is qualitatively identical to 
the case of infinite separation. For weak couplings 5^1, 
there is small but finite number of condensed counterions 



1^1 1 1 1 1 1 1 1 1 1 . 
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FIG. 2: The fraction of condensed counterions r = Zric/no 
as a function oi g = Z^Ib / ^ for different values oi 6 = . 
At low surface charge <; <C 1, the counterion distribution is 
well described by PB theory since t ^ 1. However, at high 
surface charge, correlation effects leads to a large fraction of 
counterion condensed. The condensation is first-order for 6 < 
8c and continuous for 9 > Oc- The critical point is at 6c ~ 
0.017, gc ~ 1.23, and ~ 0.43. 




attractive^ 
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g=i.4 



6=0.06 
, 6=0.02 



0.5 1 



FIG. 3; The fraction of condensed counterions as a function 
of distance. 



but the total pressure remains repulsive. For sufficiently 
high g 1, the condensation proceeds continuously for 
6 > dc and via a first-order phase transition for 9 < 9c 
at a particular value of the coupling constant go{d^9). 
We note that in this regime, the number of counterion 
condensation becomes significant. This implies that the 
mean-field repulsion is drastically reduced and the cor- 
related attractions can overcome the repulsion at finite 
distance. For 9 ^ 0.02, roughly corresponding to divalent 
counterions at room temperature, we find that the onset 
of the attraction occurs at g ~ 1.6 or surface charge of 
about one charge per 10 nm^ at a distance d = 1.5A ~ 40 
A. These numbers are order of magnitude consistent with 
computer simulations [ p^ . 

Next, we discuss counterion condensation and the to- 
tal pressure of system as a function of distance. Note 
that this scenario is more physically relevant, since sur- 
face force experiments usually vary the distance be- 
tween charged surfaces rather than changing their sur- 
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FIG. 4: The pressure profile for monovalent (Q — 0.1) and 
divalent {9 = 0.02) counterions in the case of low surface 
charges g. 



FIG. 6: First-order binding transition: for the case of trivalent 
counterion (d = 0.01), the number of condensed counterions 
(r) exhibits a discontinuous jump at a particular distance. 
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FIG. 5: The pressure profile for monovalent counterions in 
the case of moderate coupling <? ~ 1. 



FIG. 7: The presure profile for the first-order binding transi- 
tion. 



face charge densities. For low surface charges g ^ 1, 
as shovifn in Fig. ^, the counterion condensation is con- 
tinuous as a function of separation and the fraction of 
condensed counterion remains smaU but finite. We note 
that T generally increases as the distance of two surfaces 
decreases, but remains fairly constant up to A. This is 
not surprising since there is an entropy loss of the free 
counterions due to confinement. However, the pressure 
remains repulsive and shows little difference from the PB 
pressure profile, as expected. 

For sufficiently high coupling g ~ 1, we have several 
interesting regimes depending on the reduced tempera- 
ture 9 (see Fig. For 9 > 9c, the counterions condense 
continuously as the separation d decreases and the pres- 
sure of the system remains repulsive down to very short 
distances, as shown in Fig. where we have plotted the 
pressure profile for monovalent counterions {9 = 0.1) for 
different values of g. Note also that there is still a repul- 
sive barrier, which decreases with increasing g, while the 
range of the attraction is shifted to larger separations. 
For g = 1.2, corresponding to a surface charge density 
of about one charge per E ~ 300 A^, the total pressure 
becomes attractive at about d ^ 10 A. It should be noted 



that in real experimental settings, other strong repulsive 
force, such as hard-core or hydration force, that we have 
not taken into account in our model, may become im- 
portant and may overwhelm this correlated attraction 
at length scale less than ^ 20 A |jl|. This may explain 
why attraction is difhcult to observe experimentally for 
monovalent counterions. Moreover, the pressure profile 
for large separations is similar to that of the PB theory, 
except with a renormalized or effective surface charged 
density. Indeed, it is known experimentally that in order 
to fit experimental data to the PB theory, it is necessary 
to use an effective surface charge, which is always lower 
than the actual surface charge density ||l|. Therefore, 
this counterion condensation picture provides a possible 
scenario in which this phenomenon can be accounted for 
theoretically, without invoking charge regulation mecha- 
nism. 

However, for 9 < 9c the behavior of the order parame- 
ter T and the total pressure of the system are qualitatively 
different (see Fig. ^). We find that there is a range in 
the coupling constant: goo{9)/2 < g < goo{9), in which 
the order parameter displays a finite jump at a particu- 
lar separation do{g,9), and the counterion condensation 
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: 0=0.1 




FIG. 8: The pressure profile for high surface charges. Note 
that the distance at which the pressure turns attractive can 
be large compared to the Gouy-Chapmann length A. 



Finally, For g > goo{d) and 9 < 9c, the condensation 
again becomes continuous. This is because the first-order 
phase transition has already occurred at infinite separa- 
tion, in which r ^ 1. In this regime, the length scale 
at which the attraction starts to overcome the repulsion 
can be quite large (see Fig. ^). In the case of diva- 
lent counterions, we find that the onset of attraction oc- 
curs at d ~ 100 A. Clearly, the higher the surface charge 
or g, we have longer the range of the attraction; there- 
fore, together with the mechanism of fluctuation-driven 
counterion condensation, the correlated attraction may 
explain the long-ranged attractions observed experimen- 
tally. Moreover, we note that there is a qualitative change 
in the shape of the pressure: the repulsive barrier disap- 
pears. This may mark an onset of the aggregation and 
has important experimental manifestations on the phase 
behavior of the macroions. 



is first-order as a function of the separation d. Here, 
900(9) denotes the coupling constant at which the first- 
order counterion occurs at infinite separation, i.e. an 
isolated charged plate (see Sec. || and Ref. Q). Thus, 
in the limit g —^ goc{9), we must have do{g,9) 00, 
since the system is composed of two isolated charged 
plates. On the other hand, we have do{g,9) ^ as 
9 ~^ 9<x}{9)/2, because this limit corresponds to a sin- 
gle charged plate with twice of the surface charge den- 
sity, i.e. a = 2eno. This striking behavior of the order 
parameter has interesting implications for the interaction 
for the system. Indeed, for sufficiently short distances, we 
find that the first-order counterion condensation sponta- 
neously can take the system from the repulsive to the 
attractive regime, resulting in a novel first-order bind- 
ing transition. This is illustrated in Fig. ^ for the case 
of trivalent counterions at room temperature 9 = 0.01 
at 9 = 1.3, corresponding to a surface density of one 
charge per S ^ 70 nm^. (For trivalent counterions, the 
first-order counterion condensation occurs in the range of 
0.9 < g < 1.8 The corresponding pressure profile is 

plotted in Fig. ]^ which shows that the binding transition 
occurs at about d ~ 10 A. We note that an interesting 
consequence of this first-order binding transition is the 
existence of metastable states, which may have impor- 
tant manifestations in surface force experiments. It is 
easy to imagine that the system can be trapped in differ- 
ent metastable states, and therefore, hystersis may occur 
as the two surfaces are pushing in and pulling out again. 
Indeed, there are some experimental support for this be- 
havior for multivalent counterions in similar systems . 
It is important to emphasize that this interesting behav- 
ior is not included in the mean-field PB theory. Note 
also that this first-order binding transition can only take 
place at short distances. This is because do{9,9) gener- 
ally increases with increasing g, and eventually when g is 
near goo{6), the condensation occurs within the repulsive 
regime and the binding transition becomes continuous. 
Thus, direct experimental observation of the first-order 
binding transition may prove subtle. 



In summary, by incorporating the condensation driven 
by fluctuations, we show that the net pressure between 
two similarly charged surfaces becomes negative, hence 
attractions, at a length scale much longer than the Gouy- 
Chapmann length. We also predict several distinct be- 
haviors of the system, depending on the valence of the 
counterions, that deviates significantly from the classical 
theory of the double-layer interactions. While our calcu- 
lation is based on the Gaussian fluctuation theory which 
may break down for very high surface charge density, a 
complementary treatment is considered by Shklovskii |^ 
in this regime, where the condensed counterions are as- 
sumed to form a 2D Strongly Correlated Liquid. That 
theory predicts a strongly reduced surface charge and ex- 
ponentially large renormalized Gouy-Chapmann length, 
qualitatively similar to our results that for high surface 
charge most of the counterions are condensed. Moreover, 
it was shown in Ref. ||2^ that by perturbing around the 
low temperature Wigner crystal ground state, the long- 
ranged attraction persists to be operative, independent 
of the ground state. Thus, at large distances, we believe 
that our picture should capture the interaction of two 
similarly charged surfaces in the regime between where 
PB theory is appropriate (low surface charge) and the 
strong coupling limit |^, ^ . 

However, there remain fundamental issues to be ad- 
dressed in the future. For example, in real systems, the 
charged surfaces are often characterized by discrete sur- 
face charge distribution. In recent studies |33|] , it is shown 
that the counterion distribution is strongly modified if 
discreteness is taken into account. In particular, the 
counterions tend to be more "localized" near the charged 
surface. It remains to be seen how this affects the con- 
densation picture presented in this paper; it is possible 
that this effect may smooth out the first-order transi- 
tion. However, we believe that a rapid variation of the 
condensation reflecting the flrst-order transition should 
remain. 
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APPENDIX A: DERIVATION OF THE 
FLUCTUATION PRESSURE 

In this appendix, we present a detail derivation of the 
pressure arising from fluctuations. The derivative of the 
fluctuation free energy with respect to the distance d is 
given in Eq. ( p7| ) 



dd 



— J2 5iz±d/2)+ At2e-^(-) 



(Al) 



where K^e-vC^) ^ 2a'^ sec^{az) e{z) + 2 {\z\ -d/2 + ^y^ e{z) as defined in Eq. ^,C = \r/{1 + P), and Ga^lx, x') 
is the Green's function defined as the inverse operator of K-^d and satisfies the equation Eq. (p8|), which in Fourier 
space can be written as: 



Xi 



■^S{z±d/2) + K^e- 



ip{z) 



Gzd{z, z'\ q) = £b S{z - z'). 



The Green's function can be solved by standard technique [|^; first, we note that the homogeneous solutions are 
given by 



h±{z;q) = e 



1 ± — tan(Q!z) 
9 



for \z\ < d/2 and 



h±{z;q) 



It 



q {\z\-d/2 + i) 



(A2) 



(A3) 



for \z\ > d/2. We have two cases to consider: \z'\ < d/2 and \z'\ > d/2. In the former case, we split the space into 
four regions: z < —d/2, —d/2 < z < z' , z' < z < d/2, and z > d/2 and write 



G-iz,z';q) = A{z')h>{z;q), 

G^{z,z';q) = B{z')h<{z;q) + C{z')h<{z;q), 

G>{z,z';q) = D{z')h<{z;q) + E{z')h<{z;q), 

G+{z,z';q) = F{z')hl{z;q), 



for z < -d/2, 
for -d/2 < z < z', 
for z' < z < d/2, 
for z > -d/2. 



The coefficients A{z') . . . F{z') are determined by the following boundary conditions 
GU±d/2,z';q)=G^d{±d/2,z';q), 

dzG3d{z,z';q)\^^^^^2- d^G^diz, z' ■,q)\^^^j^i^ = {2/ Xd) G3di±d/2, z'; q), 

G3d{z',z';q) = G^diz' ,z';q), 
d,G3diz,z';q)\^^^, - d.Gsdiz, z'; q)\^^^, = ig- 

After some algebra, we obtain specifically 

^ h<{z')+M{q)h<{z') 
' 2q {l + a-'/q^)[l- M^{q)y 
D{z') = M{q)E{z'), 

e-grf [ (1 + 62)2 + _ 52 _ + ^2 + ] 



M{q) 



(1 -I- 62)2(1 + gAfl) + (7 + qXn){\ + 62 + q\j,){l - 62 + qXn) ' 



(A4) 
(A5) 
(A6) 
(A7) 



(A8) 
(A9) 
(AlO) 
(All) 



(A12) 
(A13) 
(A14) 
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where 7 = Xr/Xd = 2t/(1 — r). Therefore, the Green's function G'3d(x,x) for \z\ < d/2 is exphcitly given by 



(27r)2 2q 1 ^ M^q) q^ l + a^q^ 1 - M^q) 



+ 



2M{q) 
1-MHq) 



[1 + a^/q^ tan^(az)] cosh(2qz) + 2a/q tan(az) sinh(2(jz) 
1 + Q;2/g2 



(A15) 
(A16) 



Note that the Green's function is symmetric with respect to z, as expected from the symmetry of the problem. Similar 
calculation can be done for the case \z\ > d/2 and the result is 



Gr,(x,x) = 



(27r)2 2q 



^ < 1- 



1 - C{q) 6-29(1-1-^/2) [l + q [\z\ - d/2 + ^) ] 



where C{q) is given by 



g2(|z|-d/2 + e)2 
e-'J''g{d/2) h>{d/2) e-i'^ 



[h>{d/2)r 



h>{d/2) 



and 



Q{d/2) = 



[h<{d/2) + M{q)h<{d/2)][h<{d/2) + M{q)h<{d/2)] 2q 



^-^G^d[d/2,d/2-q]. 



(l + a2/q2)[i__^2(q)] 
Returning to the expression in Eq.( |Al| ), we note that it can be separated into three parts 



1 5/3AJ^A 
^ dd 

1 dp/^Ts 
A dd 

1 d(3ATc 
A dd 



IbX 



B^D 



d^Cj 

(2^ 



dz Cadiz, z;q) — S{z- d/2), 



1 f d^q f°° d 
7^7(2^ J ^^'^3d(z,z;g)0(z)— [2a^sec2(az)] 



1 



d^q 

(2^)2 



dz G3d{z,z;q)Q{z) — 



{z-d/2 + iY 



(A17) 

(A18) 
(A19) 

(A20) 
(A21) 

(A22) 



where we have used the fact that the integrand is symmetric with respect to z. Note also that there should also 
be two terms containing dQ{z)/dd and dQ{z)/dd in OfSATs/dd and df3A!Fc/dd, respectively; however, they cancel 
identically when they are added together. 

Let us first discuss Eq. ( |A20| ) ; using the identity: ^ S{z — d/2) = — S{z — d/2), and integrating by part, it can 
be transformed into 



1 dpAJ^A 
A dd 



1 



d^q dG3diz,z;q) 



IbXd J (27r) 



dz 



(A23) 



z=d/2+ 



Using the boundary condition in Eq. (A9) 

dzG3.d{z,z;q)\^^^/^+ = 9zG3rf(z, z;9)|^^^/2- + (2/Al,) G3d(d/2, d/2; q), 



and the explicit expression for the Green's function given in Eq. (A16), we obtain after some algebra 



1 dpATA 
A dd 



1 



d^q £1 



-MHq) 



qXoa + b^ib^ + qXR + qXD) 



£bXd J (2^)2 2qXD I 1 - MHq) (1 + 52)2 + _p_ + 52 + ^Xr) 

qXp b^l + b^)^ 

1 - MHq) (1 + 52)2(1 + qXj,) + (7 + gAfi)(l + 52 + qXR){l - b^ + qXn) 
d^q qM'^iq) 
(2^ l-M^q)' 



The next term, Eq. (A21), can be shown to be 
1 dl3ATB 8 (1 + 52)5 



A dd 



d^q 

£bX% 2+^(1 + 52)7 (2^ 



ad/2 



dx Gsdix, x; q) sec^ x (1 + x tana;) . 



(A24) 



(A25) 
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In evaluating the x-integral , we note that there is a nontrivial integral which involves the last term inside the bracket 
of Green's function in Eq. ( A16 ); it reads 

Q= da; sec^ a;(l + .Ttanx) { [l + (2/fc)^ tan^ x] cosh fcx + 4/fc tana; sinhfcx} , 
Jo 

where k = 2q/a and d = ad/2. Note that none of these integrals can be expressed in terms of elementary functions, 
but integrating by parts several times, one can show that the integral Q can be expressed in closed form with the 
help of the relation: 2bta.n{ad/2) = 1 — 6^, by 



Q 



16 6(gA«)2 

(1-62)25 



(1 + 5^) 



coshlqd) 



6(1 



r 



4(gA;?j 



6(1 + 62) 
2(<zAfl) 



cosh{qd) 



smh(qd). 



bcosh{qd) (l-62)(l + 62) 



'biqXR) 



d 



Xr 

:i+62 



(1 + 62) 



coah.(qd) 



smh{qd) 



Substituting this result back into Eq. (A25) and rearranging terms, we obtain 



1 d^AJ^B 
A dd 



(1 



/ 



rf2q 



£bXI 

gW2) 



d2q 



(27r)2 2q 



Q{d/2) 



J{d/2) 



(27r)2 2q 



J{dl2) 



2 

2(l-62)^(d/2) 



1 



462 _^ ^2) 



IbXI 2 +£(1 + 62) 
2 2M{q) smh{qd) 



where J{d/2) is defined by the expression 

J{d/2) 



(gAfl,)2 (1 + a2/q2) q\j, [I + / q2)[l - M^{q)] 
l+7W2(q) 2M{q) cosh qd 



1-M2(q) 1-M2(^) 

Finally, we turn to the last term in Eq. ( |Al|) , Eq. ( [A22| ). With the help of the integral 



dz ■ 



G3d{z,z;q) (1 + 62)2 



d/2 



{z-d/2 + ^y 2q 2A 



■5id/2) 



(A26) 



(A27) 



Eq. ( A22 ) can be written as 

1 dpAJ'c _ (1 + 62)2 f d2q 



A dd 



(27r)2 2q 



g{d/2) + 1 - C{q) 



462 (1 + 62) 



2+3^^(1 + 62)7 (2^)2 2q 



d2q 



g{d/2) + 1 - C{q) 



which can be combined with the expression for d(3A!FB /dd above (note that Q{d/2) cancels nicely) to yield 



IDPATb+c (1 + fo') 



,2^2 



462 (1 



A dd 

where we have defined the following dimensionless integrals 

rf2q AtiXb. \ 2{l-b^)J{d/2) 



£bXI AttXr " ' £bA2j 2+3^(l + 62)4^Ai 



(l2 



I2 = 



(27r)2 2q 
d2q 47rAi 



2M{q) sinh(gd) 



_ (gAfl)2 (1 + a2/^2) (1 + a^/q') [1 - M^{q)] 

l^lL{q)~J{d/2) 



(27r)2 2q 

With some straightforward but tedious algebra, they can be cast into more explicit form: 

2x [(1 - 62 + a;)2 + 7(1 + 62 + x){l - b^ + x) - b^ { 462 _^ ^2 



(A28) 

(A29) 

(A30) 
(A31) 



I2[d/Xi 



dx ■ 



dx 



(462 + 2.2) [1 _ M^ix)] [(1 + 62)2(1 + a;) + (7 + a;)(l + 62 + x)(l - 62 + x)] 
2xMHx) {{1 - b^ ^ x) [(l-62+x)+7(l + 62+a;)] - ( 62 + x) (462 + x2)} 
(462 + x2) [1 - M^{x)] [(1 + 62)2 + - 62 - x)(l + 62 + x)] ' 



(A32) 
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and 



/>CXD 

l^d/XR] = - / dx- 
Jo 



dx 



2jxb'^ 

M^{x)] [(1 + 62)2(1 + x) + (j + x){l + b'^ + x){l - + x)] 
2xM'^{x) [x + 7(fe2 +a;)] 

[1 - M^{x)] [(1 + 62)2 ^(1 _ 52 _ + ^2 _^ 3.)] ' 



(A33) 



where we have made a change of the i ntegr ation variable x = qXr. Now, observe that the first term in Eq. ( A2E ) 
cancels precisely the first term in Eq. ( A24 ). Therefore, combining the two expressions, we obtain Eq. ( |2g| ) for the 
fluctuation pressure. Similarly, 2i[d/\ft defined in Eq. ( ^6|) can be expressed as: 



Ji [d/\ 



d^q AttXr 



qXR, [(1 + 62)2 +27(1- 62 +gAfl,)] 



(27r)2 2q (l + 62)2(l + gAfl) + (7 + <zA«)(l + 62 + 9Afi)(l 
d2q 47rAfl^ qXRM^iq) ( (1 + 62)2 ^ 37(1 - 62 + gA^) 



62 + -7Ak) 



(27r)2 2q 1 



A12(q) t (1 + 62)2(1 + qXj,) + + qXj,){l + b^ + qXj,)(l 
_ (1 + 62)2 ^ 2qXR{l - 62) + 2(gAfl.)2 + 27(1 - 62 - gAj,) 
(1 + 62)2 + ^(1 _ 52 _ g^^j(i + 52 + qXR) 



62 + gA;^) 



(A34) 



Note that the first term in this expression is logarithmi- 
cally divergent. 
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